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Extensive Monte Carlo results are presented for the structure of a bottle-brush polymer under 
good solvent or theta solvent conditions. Varying the side chain length, backbone length, and the 
grafting density for a rigid straight backbone, both radial density profiles of monomers and side 
chain ends are obtained, as well as structure factors describing the scattering from a single side 
chain and from the total bottle-brush polymer. To describe the structure in the interior of a very 
long bottle-brush, a periodic boundary condition in the direction along the backbone is used, and 
to describe effects due to the finiteness of the backbone length, a second set of simulations with 
free ends of the backbone is performed. In the latter case, the inhomogeneity of the structure in 
the direction along the backbone is carefully investigated. We use these results to test various 
phenomenological models that have been proposed to interpret experimental scattering data for 
bottle-brush macromolecules. These models aim to extract information on the radial density profile 
of a bottle-brush from the total scattering via suitable convolution approximations. Limitations 
of this approach and the optimal way to perform the analysis of the scattering data within this 
approach are discussed. 





I. INTRODUCTION 



Recently __there_has__been a great experimental (see 
and theoretical 

1^4 liij III IM Pi 11^ 123 

[3ll . [32 . [33l . [SJ, [3^ interest in the 
conformation of so-called bottle-brush polymers. Such 
polymers consist of a long (flexible) main chain, at which 
many flexible (shorter) side chains are densely grafted, 
such that an overall shape of a worm-like cylindrical 
brush results P, i, i, 11 Q- Synthesizing such poly- 
mers with suitable characteristics, materials can be pre- 
pared whose properties can be adjusted by external stim- 
uli, such as the solvent quality, pH, or temperature, and 
this fact makes such bottle-brush polymers interesting 
for various applications [s^ [s^l- For controlling the 
properties of such bottle-brush polymers, a good theo- 
retical understanding of their structure and conforma- 
tion, as a function of control parameters such as the 
chain lengths of the main chain and the side chains, 
their grafting density, and the solvent quality, is manda- 
tory However, despite the longstanding theoretical ef- 
fortsj^, 10, n, 12, 13, 14, 15, 16, 17, 18, 19, M [H, 

111, 111, ll 111 IM |27l, 12a, l2g, |3d, IM ^ this 
problem is still incompletely understood. While one has 
various scaling predictions (see 3J| for a recent review) 
and treatments based on the strong stretching limit of 
self-consistent field theory exist since a long time (e.g. 
Refs. [sl. [9I [1 ^ ITl | ). recent simulations [s^l and also the 
experiments [3|, \% |5| indicate that the regime where these 
theories become accurate would require side chain lengths 
that arc hardly accessible either by experiment or by sim- 
ulation. As a consequence, the theoretical guidance for 
the interpretation of extensive experiments by combined 
light and small-angle neutron scattering analysis 0, 0, Q 
is still incomplete. 



In the present work, we make a contribution to clar- 
ify this problem by extensive Monte Carlo simulations 
of bottle-brush polymers [s^, [11] using the "Pruned- 
Enriched Rosenbluth Method" (PERM algorithm) [Hi 
to obtain the relevant information on the conformation 
of bottle-brush polymers under various conditions, that 
are needed to test the phenomenological models used to 
interpret the experimental scattering data As 
was shown in [3J| the PERM algorithm is very power- 
ful to obtain a wealth of simulation data for the case 
of side chains grafted to a hard rod (a generalization 
of the algorithm to a flexible backbone is far from triv- 
ial), representing a strictly rigid backbone polymer. This 
idealization describes a real bottle-brush chain only lo- 
cally. However, all theoretical models used for the anal- 
ysis of experiments [1, [3, [E] determining the structure of 
a bottle-brush do contain the rigid backbone as a spe- 
cial case. It is this case for which we can undertake a 
stringent test of physical model assumptions underlying 
the analysis of experimental data. Of course, there is 
no reason to assume that a model that already fails for 
the (simpler) rigid backbone case should be accurate for 
bottle-brushes with flexible backbones. 

In Sec. II, we introduce our model and recall the most 
basic facts on our simulation method and define the quan- 
tities that are studied. In Sec. HI, we give a comprehen- 
sive overview of our results on various physical properties 
of the bottle-brush polymers. Sec. IV then is devoted 
to the problem relevant for the interpretation of the ex- 
periments, namely the test of theoretical models used 
in [1, [3, [1, [^ for our system: Note that unlike the ex- 
periment, we can extract information of radial density 
profiles and geometrical characteristics of individual side 
chains directly from the simulation, simultaneously with 
but independent of the information gathered on the scat- 
tering functions, and thus a stringent test of the proposed 
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relations between these quantities is possible. Sec. V then 
summarizes our conclusions. 



II. MODEL AND SIMULATION 
METHODOLOGY 

As in Ref . [13] , we use as a coarse-grained simple model 
of flexible polymers in solution, the self-avoiding walk on 
the simple cubic lattice. Each lattice site can be occupied 
by a single effective monomeric unit only, and this ex- 
cluded volume interaction corresponds to polymer chain 
conformations under good solvent conditions [3^ H^, [4l[ . 
Introducing an energy parameter e that is won if two ef- 
fective monomers occupy nearest neighbor sites, one can 
describe variable solvent quality in this model simply by 
varying the temperature T: the Theta temperature O 
where this attraction approximately cancels the excluded 
volume repulsion, in the sense that the mean square gy- 
ration radius {R'^)t of a chain scales linearly with the 
chain length N, apart from logarithmic corrections [s^ . 
occurs for ge = exp(-e/A:i3e) 1.3087 [38]. We shall 
present results both for T = 8 and for T — > oo {where 
q = exp(— e/fcsT) — 1, and hence only the excluded vol- 
ume interaction is present}, in view of the fact that most 
cases of experimental interest will be somewhere in be- 
tween these limits. 

Following Ref. [13] we take the rigid backbone along 
the z-axis of our coordinate system. Using even values 
for the length Lb of the backbone, measuring all lengths 
in units of the lattice spacings, grafting sites z = 1, 
z = Lb at the end of the backbone are labeled as si, 
sites adjacent to the ends {z = 2, z = Lb-i) as S2, and so 
on, until the center of the bottle-brush, sites z = Lb/I, 
z = Lb/2-\- 1 being denoted as Si(,/2- Of course, a depen- 
dence of properties of a side chain on the coordinate 
occurs only when we consider bottle-brush polymers with 
free ends in the z-direction, while no dependence on Sk 
occurs if we choose periodic boundary conditions (pbcs) 
in z-direction such that the grafting site z = Lf, is nearest 
neighbor of z = 1: in this case full translational invari- 
ance in z-direction holds, and the distribution function 
of the monomers Pz' (z — z') of a side chain grafted at z' 
must be symmetric around z', Pz'iz — z') = Pz'iz' — z). 
This symmetry property does not hold only for the distri- 
bution function of all the monomers that belong to that 
side chain, but also for individual monomers i = 1, . . . , N 
along the side chain, in particular for the chain ends. Also 
the average z-coordinate of the center of mass of the side 
chain coincides with z'. None of these symmetries carry 
over to the case with free ends, of course; in the latter case 
the monomers of the side chain can have z-coordinates 
in the region —N + l<z<Lb + N (the boundaries of 
this interval do not occur in practice, of course, it would 
require that a side chain grafted at z = 1 or z = Lf, is 
linearly stretched out in the —z or +z-direction, respec- 
tively) . 

When one considers properties of individual side 



chains, which are stretched away from the backbone, 
two non-equivalent directions x, y need to be distin- 
guished [34| : defining the vector toward to the center 
of mass (CM.) of a chain from its grafting point z' as 
{X z' ,Yzi T Z z') in a fixed laboratory frame, for a partic- 
ular configuration of the side chain, we can define the 
X-axis along the vector (X^/,Yz'), and require that the 
j/-axis is perpendicular to the x-axis and also lies in 
the Xz' — Yz> plane. Since for a densely grafted bottle- 
brush polymer strong stretching of the side chains is ex- 
pected d, i, O, [H [13, [H, [11 this distinction allows 
to compute linear dimensions of the side chains in the 
direction along which the chain is stretched, and perpen- 
dicular to it. 

In practice, side chain lengths up to TV = 50 were 
considered, while choices Lb = 32, 64, 128 and 256 as 
well as two values of the grafting density, ct = 1/2 and 
tj = 1, were considered. A distinctive feature of our 
implementation of the PERM algorithm is [13] that in 
one run one gets information on properties for all in- 
teger values of N from = 1 up to A'max {which in 
our case was chosen to be Mnax = 50, so the largest 
polymer simulated had a total number of monomers 
A^tot = Lb + LbaN = 256 + (256/2)50 = 6656, since for 
Lb = 256 the case ct = 1 no longer was feasible}. For de- 
tails on the implementation of the PERM algorithm for 
bottle-brush polymers, we refer the reader to Ref. [34] , 



III. STRUCTURAL PROPERTIES OF BOTTLE 
BRUSHES: THE EFFECT OF CHAIN ENDS 

In this section we will look at the difference in structure 
at the free ends of the backbone, where we can expect 
to find star-like conformations for the side chains, and 
the central part of the backbone which will be brush- 
like. Comparing conformations for bottle-brushes with 
free ends to those where pbcs are employed along the 
rigid backbone, we can also find out to what extend the 
free ends influence the average structure of the brush and 
its side chains. We will perform this comparison for good 
solvent as well as for theta-solvent conditions. 

Fig. [1] presents our data for the perpendicular part of 
the mean square gyration radius, (R^ where j_ = 

Rg ,j. + Rg y, aivl X 8.11(1 J/- c o mp oucut s refer to "mea- 
surements" taken in the laboratory system with fixed ori- 
entations of the coordinate axes along the axes of the 
simple cubic lattice. One sees that (i?^ ^) for the graft- 
ing density cr = 1 is always larger than for a — 1/2, while 
the dependence on backbone length is almost invisible. 
In the good solvent case, data for the decade 5 < A^ < 50 
are compatible with a power law increase, but the expo- 
nent is far too small in comparison with the prediction of 
the scaling theory {for large enough A^ and high grafting 
density one expects 0] {R^g J/N^" cx N^Hi-'^)/ii+-^) ~ 
^0.305 Tff]n[iQ ^i^g effective exponents that one can read 
off from Fig. [T^ are only about half of this value}. In- 
terestingly, also in the 0-solvent case one observes an 
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FIG. 1: Log-log plot of the rescaled mean square gyration ra- 
dius perpendicular to the backbone, (_Rg x)/-^^" of the whole 
bottle-brush versus the side chain length, for a good solvent 
where u = 0.588 (a) and a 0-solvent where v = 0.5 was taken 
(b). Two choices of a and four choices of Li, are included, as 
indicated. All data are for bottle-brushes with free ends. 



increase of (i?^ j_)/N'^'^ = (i?^ j_)/N with increasing side 
chain length TV, but there clearly occurs curvature on 
the log-log plot, and thus already the data indicate that 
the asymptotic region where power laws and scaling con- 
cepts apply is not reached. Analogous data have also 
been taken for the model with pbcs, but the data are al- 
most indistinguishable from the free end case, and hence 
not shown here. 

Fig.[2]now turns to the linear dimensions of side chains, 
using a coordinate system where the x-direction is de- 
fined from the direction of the vector through the back- 
bone and the CM. of each side chain, and perpendicu- 
lar to the backbone direction in each configuration (see 
Sec. II), and also different grafting sites are distinguished, 
for a relatively short backbone length. Lb = 32. As ex- 
pected (Fig. ^^), the stretching of chains grafted near 
the free ends (si) in radial direction is weakest, because 
they acquire a noticeable component in the z-direction 
(Fig. [2j;) . These effects rather quickly get weaker when 
the grafting site is farther away from the chain ends, 
and even for a short backbone {Li, = 32) the chains 
near the center almost behave like chains in the bulk 



of a very long chain (which is modeled by eliminating 
end effects through the choice of pbcs). We do not 
have such an obvious interpretation for the weak (but 
for the backbone ends clearly non-monotonic) variation 
of {RlJ/N^", however. 

It is interesting to contrast these results to the 0- 
solvent (Fig. [3]) . In this case the inhomogeneity caused 
by the presence of free ends of the backbone is much 
weaker, the differences with respect to the p.b.c. case 
are much less significant. However, a rather strong ef- 
fect of the inhomogeneity in the z-direction is seen when 
one considers the corresponding components of the mean 
square backbone to end distance of the side chains, and 
this effect is present both in the excluded volume case 
and in the 0-solvent case (Fig. 2]). 

The next question we ask is the following: how likely 
is it that monomers (or chain ends) are not in the re- 
gion 1 < z < Lb where the grafting sites are? Fig. [S] 
shows also from this criterion that in the good solvent 
case the bottle-brush is more extended in the z-direction 
than in the 0-solvent case. Even for short backbones 
{Lb ~ 32) for 0-solvents bulk behavior is reached, while 
for good solvents there is still some systematic depression 
in the center (z = Lb/2 = 16). We note, however, that 
for larger Lb such as Lb = 64 (to save space these data 
are not shown) bulk behavior is reached for a significant 
range of z in the center of the bottle-brush. 

Another quantity that shows that side chains near the 
backbone ends tend to orient much more along the z-axis 
in the good solvent case rather than in the 0-solvent case 
is the distribution P{9) of the angle between the vectors 
towards the center of mass of each side chain and the 
z-direction (Fig. [6]). One should note that angles 9 near 
6 = tt/2 characteristic for chains stretched away from the 
backbone in perpendicular direction, dominate only in 
the center of the backbone, while angles near 6 — it /A and 
37r/4 make a substantial contribution near the backbone 
ends. For the considered side chain length, this effect 
dies out after a few monomeric distances away from the 
backbone ends, however. For 0-conditions (Fig. [6)3) this 
behavior is only found close to the chain end. Chains 
grafted already five monomers away from the backbone 
end show no tilting like for the good solvent case. The 
average angle remains at 7r/2, but the distribution gets 
broader and asymmetric with a heavy tail towards the 
adjacent chain end. 

The data shown in Figs. [2][6] are readily accessible in 
simulation, but not easy to access experimentally. They 
help, nevertheless, to develop a complete picture of the 
structure of bottle-brush polymers and clarify the side 
chain conformations. Quantities, that experimentalists 
try to extract from their studies are accessible to the sim- 
ulations as well, of course. Such quantities are the radial 
distribution p{r) of the monomers and pdr) of the chain 
ends (shown in Figs. [3 [8] for Lb = 64). Note that due to 
the discreteness of the lattice, the number of monomers, 
N{r) and Ne{r), in the interval [r,r + dr] are normal- 
ized i.e. p{r) = N{r)/Nr and pdr) — Ne{r)/Nr where 
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FIG. 2: Log-log plot of the rescaled mean square gyration radii of the side chains, {R%^x) /N'^" (a), {F!?gc,y) / N'^" (b) and 
{Rgc,z)IN'^^ (c) versus the side chain length N, for the good solvent case, Lj = 32, cr = 1 and various choices of the grafting 
sites, as shown by the coordinate Sfc (cf. Sec. II for explanations). The full curves show analogous data for the case of pbcs. 



Nr is the number of lattice points with a distance to the 
backbone lying in the interval [r, r + dr] . For comparing 
data of different chain lengths, normalization conditions 
^A^(r) — N and X)-^e(^) ~ have been imposed. 

r r 

Similar data have also been generated for Lt = 32, but 
the differences to those shown are only small, and there- 
fore need not be discussed here. Figs. [7] and [S] reveal 
that neither p(r) nor Pe{r) are sensitive to the effects 
of the free ends: for the chosen i^, much longer side 
chains would be required in order that effects due to 
the crossover from bottle-brush to star polymer behav- 



ior come into play. While in the good solvent case the 
chain ends are typically farther away from the backbone 
than in the 0-solvent case, the qualitative behavior of 
p{r) and pe (r) does not depend on solvent quality much. 
Furthermore, it is gratifying to note that these data are 
qualitatively rather similar to the corresponding Molecu- 
lar Dynamics results of Murat and Grest [3] for a bead- 
spring off lattice model of flexible side chains tethered to 
a straight line. This similarity reinforces our view that on 
a coarse grained level, the present lattice model should 
yield useful results. 
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FIG. 3: Same as Fig. [S] but for the case of the 0-solvent. 



IV. SCATTERING FUNCTIONS FOR 
BOTTLE-BRUSH POLYMERS AND THEIR 
THEORETICAL MODELING 



total scattering function S'w('z) for the bottle-brush poly- 
mers, both for good solvent and 0-solvent conditions. 
Here Sv^iq) is defined as 



Let us now turn to a discussion of experimentally ob- 
servable information on the structure of a bottle brush 
polymer. In experiments one has to infer the structure 
from scattering data [1, 0, Q employing suitable model 
assumptions on the structure. In the simulation we ob- 
tain both, the scattering data and the underlying struc- 
tural properties described in the last section indepen- 
dently, and therefore are able to test theoretical models 
suggested to link the two. Fig. [S] presents our data for the 



SA,) = ^f:f:icirMr,)f-^^^ (1) 
Mot ^ ^ q\ri - rj\ 

where c{fi) is an occupation variable, c{ri) = 1 if the site 
fi is occupied by a bead, and zero otherwise. Note that 
an angular average over the direction of the scattering 
vector q has been performed, and the sums run over all 
monomers (all side chains and the backbone). 

Surprisingly, our data are qualitatively very similar 
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(a) (b) 




FIG. 4: Log-log plot of the mean square backbone-to-end distance of the side chains versus the side chain length A'', for Lb = 32 
and a = 1. Panels (a,b) refer to the good solvent, panels (c,d) to the 0-solvent case. The components shown are {R^)/N'^'^ 
(a,c) and {Ri)/N^'^ (b,d). Symbols denote different coordinates Sk along the backbone, while the full curves show the analogous 
result for pbcs. 



to the corresponding experimental data (see e. g. Fig. [4] 
of 0)j although the latter refer to a polymer with a flex- 
ible backbone, unlike our simulations. As always, the 
limit g ^ of Svi{q) reflects the total number A/tot of 
scattering monomers, and the leading deviation from it 
is described by the total gyration radius, 

^wCg) «Mot[i-g'(i?g>/3] . (2) 

This behavior is shown by the flne-dotted lines in Fig. [5] 
for the case of grafting density a = 1. The q range over 
which this approximation agrees with the scattering data 
increases with increasing ratio of side chain length to 
backbone length, N/Lf,. Of course, more interesting is 
the behavior at larger q, where Eq. ([2]) is no longer valid. 
The region where {q) is strongly curved and decreases 
rapidly (0.1 < g < 0.3, in our case) has contributions 
from the conformation along the backbone (rigid rod in 
our case which should show up as a behavior Sw{q) ~ g"^ 
for longer backbones) and from the scattering from the 
cross section through the cylindrical bottle-brush, and 
needs to be related to data such as shown in Figs. [3 [S] 



The q range near q — 1 reflects the self-avoiding walk 
structure q~^^'^ before it is affected by the local pack- 
ing of monomers on the lattice at still larger q, and in 
real systems reflects local properties such as the persis- 
tence length of the flexible side chains, possible scattering 
from side groups, etc. This non-universal regime hence is 
less interesting. From this discussion of the total struc- 
ture factor we can already conclude that it is the g-range 
0.04 < q < 0.5 which for our model contains the impor- 
tant information about the structure of the brush. 

One advantage of our simulations is that we can ob- 
tain scattering contributions from different parts sepa- 
rately. E. g., we can isolate the scattering from the back- 
bone fFig.fTUb) and from the scattering of the side chains 
(Fig. [TUb.c). In our case, where the backbone is a rigid 
rod where just the subsequent sites i = 1, 2, . . . , are 
taken by monomers, Eq. ([T]) becomes 



Lb Lb . / 

z=l j=l ^ 
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FIG. 5: Density distributions of all the monomers, p{z), (a,c), and of the chain ends, pe{z) (b,d), plotted vs. z for Lt = 32, 
a = 1, good solvent conditions (a,b) and B-solvent conditions (c,d). The distributions are normalized by choosing p(z) — ric 
and Pe{z) = ric, where ric is the number of side chains {ric — crLb). Four chain lengths are shown, as indicated. Note that 
in the p.b.c. case we trivially have p{z) = Pe{z) = 1, 1 < 2 < Lb, for the chosen normalization. 



Noting that the distance |j — «| = occurs Li, times, 
while the distance \j — i\ — 1 occurs 2(Lf, — 1) times, 
the distance \j — i\ — 2 occurs 2{Lb — 2) times, etc., we 
conclude that 

+ f EV.-.)5^. (4) 

The factor 2 accounts for the fact that both positive and 
negative differences k = j — i occur, and the extra — 1 
corrects for over counting in the term k = 0. 

In the hmit where Lb 00 and qLf, is of order unity, 
the sum in Eq. ^ can be transformed into an integral, 
to find 

s,(,).l r'!!Iil,,-4!^l!(|^ii« (5) 

q Jo t q Lb 

Eq. ([5]) is nothing but the well-known scattering function 
of an infinitely thin rod of length Ly, with a continuous 
mass distribution along the rod [i^, . 



According to Eq. ([T]) , the scattering function of all side 
chains is given by 

When we add Sh{q) and Ss{q) with their relative weights, 
see Figs. fTOk.b. we do not recover S'w(9) strictly, due to 
interference effects in the scattering from monomers in 
the side chain and in the backbone. Such interference 
effects normally are neglected 0, 0, S [1) . Taking the dif- 
ference Sbs = WtotS^iq) - LbSbiq) - iVneS',(9)]/(27Vtot) 
, we can test for the importance of such interference ef- 
fects as shown in Fig. [TlJ Indeed, we do find that such 
interference effects are present although only at a percent 
level. 

We now turn to a discussion of the scattering from 
the side chains, which clearly dominates the scattering 
intensity in all cases of practical interest. For scatter- 
ing wavenumber q in the range q{Rg) 3> 1 this scatter- 
ing should be dominated by the cross sectional structure 
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FIG. 6: Distribution P(6) of the angle 6 between the vectors 
towards the center of mass of each side chain and the direction 
of the backbone, for Li, = 32, cr = 1, A'^ = 50, the good solvent 
case (a) and the 0-solvent case (b). The different symbols 
indicate different positions along the backbone, as indicated. 
The corresponding distribution for the p.b.c. case agrees with 
the Cifi/ci7 curves. 



FIG. 7: Radial distribution function p(r) of all monomers (a) 
and radial distribution function Pe(r) of chain ends (b), plot- 
ted versus r for Lj, = 64, a = 1, good solvent conditions, and 
four values of the side chain length N, as indicated. Symbols 
show our results for free ends, while curves show correspond- 
ing data for the case of pbcs. 



of the bottle-brush. In the analysis of the experimental 
scattering data one has to assume that one can determine 
the cross-sectional contribution by a factorization 

S'w(g) = Shiq)S^siq) , (7) 

where S^sil) is interpreted as the cross section structure 
factor. Such decoupling approximations seem to be suc- 
cessful for worm-like micelles [i^. In the literature, Sh{q) 
is modeled by a superposition of rigid rod and worm-like 
chain form factors, needed to account for backbone bend- 
ing 6]. In our case we can take Eq. ([7]) simply as a def- 
inition of Sxsil) using Sb{q) which is known exactly for 
our case (see Eq. (jj])). 

The cross sectional scattering is then assumed to be 
obtainable from a rotationally averaged two-dimensional 
Fourier transform of the radial density distribution. 



S^s{q) = 




Here C is a normalization, and the indices T and q in- 
dicate a thermal average and an average over the unit 
circle in two dimensions. This is further approximated 
by neglecting correlations in the radial density fluctua- 
tions 

(p(r)p(f'))T = (p(r))T(p(r'))T p{r)p{r') (9) 
to obtain 



5'xs(g) = 




With the proper normalization this yields 

I Jo arrp{r)\ 

where Jo{r) is the zeroth order Bessel function of the 
first kind. With the approximations underlying Eg. (111!) 
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FIG. 8: Same as Fig,. [71 but for 0-solvent conditions. 



the experimental cross section structure factor can be 
inverted to obtain the radial density distribution 
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FIG. 9: Log-log plot of the scattering function of the whole 
bottle-brush polymer, S'w(g), in a good solvent (a) and a O- 
solvent (b) versus q. All data are for the case of free ends, 
N — 50, while data for two choices of a and three choices of 
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theoretical power laws for flexible chains and intermediate q- 
values, S^{q) oc f?"^/"' and (a) u = 0.588 or (b) u = 0.5 (b), 
respectively. Dotted curves are given by Eq. Q for a = 1. 



1 

2^ 



-S^M]'^^MQr)qdq. (12) 



In the analysis of experimental data, different plausi- 
ble assumptions for the radial density profile were used, 
guided by the assumed similarity to the scattering from 
worm-like micelles. Rathgeber et al. iS'] propose to use 
the following empirical function 



,(g)=const[^:^i%^exp(-g2,V2)]^ 

qRc 



(14) 



where Rc is an effective radius, and s is an effective width. 
This is equivalent to assuming a radial density profile 
which is a convolution of a step function and a Gaussian 



.9(0 



for r < Rc 

{1 + exp[(r - i?s)/(Ts]}"i for r > R, 



(13) Pxs(f)=PO 



Here Rc is an inner radius, up to which p{r) is a constant; 
then there is a power law decay, described by an exponent 
k, up to some outer radius Rs, then a fast decay to zero 
(over the range as) follows. The constant a is fixed by 
the condition that g{r = Rc) is continuous, so Eq. 
involves the four nontrivial fitting parameters Rc, k, Rs 
and as- Zhang et al. [3] assume a form for the cross 
section structure factor in terms of the first order Bessel 
function Ji{x), 



d^f' [1 



Rc)] exp 



(f-f')^ 



2s2 



(15) 

Here we have three free parameters, po, Rc and a. Again, 
Rc is a measure of the range over which the density profile 
is assumed to be flat in the core of the bottle-brush 

When we look at the density profiles determined di- 
rectly from simulations (see Fig.fT^). however, we recog- 
nize that there is no convex region in the interior of the 
bottle brush, even at a grafting density of one which is the 
limit of what is typically reached in experiment. There- 
fore, the comparison with worm-like micelles is mislead- 
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ing, and we suggest to use an alternative form of fitting 
function for the radial density 



h{r) 



1 + (r/ri)^i 



exp[-(r/r2)' 



(16) 



where a is the grafting density and ri and r2 are the 
length scales for the algebraic decay close to the back- 
bone and the exponential cutoff at larger distances (i.e., 
we expect ri << r2 in the course of our fit analysis), and 
xi and X2 are the corresponding exponents. Taking into 



account the predictions of scaling theory 3J| we can fix 
the first exponent xi = {Sv — l)/2v. So again we are 
using three fit parameters. In Fig. [12] we show that this 
assumed form for the radial density is able to fit the simu- 
lation data perfectly over almost six orders of magnitude 
in density for both, good solvent and 0-solvent condi- 
tions. The parameters of the shown fits are ri = 0.49, 
xi = 0.65 [v = 0.588), rs = 10.67, and X2 = 2.80 for 
the good solvent case, and ri = 1.19, xi — 0.5 [v = 0.5), 
r2 — 7.13, and X2 = 2.18 for the 9-solvent case. Within 
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FIG. 11: Sbsiq) = [NtotS^{q)^LtSb{q)~NncSs{q)]/{2Ntot), 
case (a) is for good solvent conditions, case (b) for B-solvent 
conditions. 



the range of backbone lengths studied, the radial den- 
sity profiles agree, with some statistical fluctuations vis- 
ible for the good solvent data and the longest backbone. 
Lb = 128. 

Let us now turn to a discussion of the possibility to 
deduce the radial density profile from the cross sectional 
structure factor as defined in Eq. ([7]). Figs. [T5b andfUk 
show fits to the radial density profile using the functional 
forms g{r) and h(r) defined above. The form g(r) sug- 
gested by Rathgeber el al. jSj is able to fit the radial 
density well over about 3 orders of magnitude with pa- 
rameters Rc = 0.3, k — 0.65, Rs ~ 10.5 and a-g = 1.90; 
h{r) fits over the complete range, as discussed above. 
When we Fourier transform these functions according to 
Eq. (dJ) and compare with the cross sectional structure 
factor (full lines in Fig. [T5b and I14bl we see that the 
transform only describes the scattering data well up to a 
momentum transfer value of about q ~ 0.08. This is only 
slightly larger than the range over which one only sees 
the scattering from the large-scale structure of the bottle 
brush (Eq. ((2)), which fits the cross sectional structure 
factor up to a momentum transfer value of g « 0.04, as 
shown by the dotted line in Fig. [Hb . This regime then 
is basically determined by the normalization of the ra- 
dial density distribution. Using an iterative optimization 
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FIG. 12; (a) Radial distribution function p(r) plotted versus 
r for side chain length A'' = 50, three choices of backbone 
length Lj, as indicated, and the grafting density cr = 1 for 
good solvent conditions, (b) Same as (a) but for G-solvent 
conditions. Parameters of the fit function h(r) are quoted in 
the text. 



procedure [4J| we can also find the best fit of the Fourier 
transform of the radial densities to the cross sectional 
scattering shown by the curves indicated as Qsix) (fit pa- 
rameters are Rc = 0.3, k = 0.65, Rg = 7.5, and Us = 2.8) 
and hs{r) (fit parameters are ri = 0.49, xi = 0.65, 
r2 = 8.20, and X2 = 1.80) in Figs. [TSb and Eb, where 
we extended the fit up to g « 0.4. When we then look at 
these functions in real space in Figs [T5b and [Tlk . we see 
that they are a rather poor fit to the radial density pro- 
file. The function hm{r) (fit parameters are ri — 0.49, 
xi = 0.65, r2 = 10.20, and X2 = 2.80) in Fig. El will be 
discussed later in the text. 



Using the functional form of Eq. (|T4)) we can directly 
fit the data in q-space and then transform into real space. 
This is shown in Fig. [15] for the good solvent case and in 
Fig. [m for the 8-solvent case. Looking at Figs. fTSb and 
[TBb one first has to comment on the fact that the full 
curves in both figures do not agree with the data given 
by the symbols. For these curves, the scattering data 
in the q-range [0,27r] were Bessel transformed into real 
space and back again. The overestimation of the real 
scattering for momentum transfers larger than about 0.1 
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FIG. 13: (a) Radial distribution function p(r) plotted versus 
r for side chain length TV = 50 for good solvent conditions. 
Parameters of the fit functions g{r) (best fit to p{r)) and 
gs{r) (Fourier transform of best fit to 5'xs(fi')) are quoted in 
the text, (b) The corresponding cross section structure factor 
S^s{q) ~ Sw{q)/Sh{q) plotted in the representation qS^siq), 
vs. q. The two curves correspond to the two curves in part 
(a). 



indicates that there is intensity in the modes for q- values 
larger than 2t: virhich is aliased into the studied range. 
Hovifever, looking at the direct transform of the scatter- 
ing data into real space (pis'* in Figs. [TsT a) and frBFa)) 
one sees that this is not a relevant numerical problem. 
Assuming the whole displayed q-range to be relevant for 
the determination of the radial density profile leads to 
the prediction of a highly oscillatory non-positive radial 
density. Similarly, when we try to fit the scattering data 
beyond a q- value of about 0.4 by the assumed functional 
forms g{r) and h(r) we obtain unphysical radial density 
profiles. Constraining the fit with the functional form 
of Eq. (|14p to the q-range below 0.4, however, also does 
not lead to a satisfactorily prediction of the radial den- 
sity as can be seen in Figs. [TSk and [T6k . The assumed 
convex shape of the radial density leads to an overesti- 
mation of the density in the interior of the brush and a 
compensating underestimation in the outer parts. 

Summarizing this discussion we have to conclude that 
there is only a small range of momentum transfers where 
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FIG. 14: Same as Fig. 1131 but the fit functions are h{r) and 
hs{r). Furthermore another fit is included which is extended 
up to the maximum in part (b) of this figure, h,n{r). 



the analysis using Eq. (jlip may be warranted. This range 
extends at most to the position of the maximum in the 
plot of qSxsiq) vs. q. In this q-range, one should employ 
a concave fitting function like the empirical law given 
by the function h{r) above and not the convex forms 
usually assumed for the inner part of the brushes. The 
grafting densities typically employed in experiment are 
not high enough to lead to a radial density which resem- 
bles a filled cylinder with a smeared out interface to the 
solution when one works at good solvent or 0-solvent 
conditions. This assumption may be valid working in 
poor solvent, a regime which was not accessible to us us- 
ing our simulation approach. When we perform a fit to 
the cross sectional scattering only for momentum trans- 
fers smaller than the peak position in the plots of qS^sil) 
vs. q, we obtain the function h„i{r) included in Fig. 1141 
We can see that this is a good representation of the radial 
density down to values of about p = 0.01. 

As a final result let us discuss the cross sectional radius 
of gyration of the brush defined by 



(17) 



Table 1 gives the resulting radii of gyration for the dif- 
ferent fitting functions and procedures employed and for 
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FIG. 15: (a) Monte Carlo data for the radial distribution 
function p(r) of the monomers plotted versus r for the side 
chain length N = 50, a — 1 the good solvent case, and back- 
bone length Lt = 32. The full curve labeled piV (r) shows the 
result of Bessel transforming (Eq. (|12|l ') the simulation data 
for the scattering function into real space. Pxs('") shows the 
prediction for p(r) obtained from fitting Sxa(9) using Eq. (|14p 
(fit parameters are Rc = 1.0, and s = 3.67). (b) The cross 
sectional scattering qSxs{q) is plotted vs. q. Symbols are data 
points, sis' (full line) is the Bessel transform of the full curve 
in part (a) which should ideally coincide with the symbols 
(see text). 5'xs (dashed line) is the best fit of Eq. (|14|l to the 
data. 



TABLE I: Results for the cross sectional radius of gyration 
(see text) for the different fitting procedures and both solvent 
conditions. 





h{r) hm(r) hs{r) 


g{r) gs{r) 


pUr) 


Rgc (good solvent) 
Rgc (0-solvent) 


7.83 7.49 7.64 
6.19 5.87 5.98 


7.94 7.79 
6.39 6.27 


5.24 
4.08 



the cases of good solvent and 0-solvent. All fitting pro- 
cedures reproduce the shrinking of the brush going from 
good solvent to 0-solvent condition. The results using 
the functional forms h{r) and g{r) agree well with each 
other and also the suggested fit analysis of the scatter- 
ing yielding function h„i{r) results in only 4% deviation 
from the true value. The fits using Eq. p4p . however, 
underestimate Rgc by about 33%. 
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FIG. 16: Same as Fig. [15] for 0-solvent conditions, (fit pa- 
rameters are Rc = 1.0, and s — 2.88) 



V. CONCLUDING REMARKS 

In this paper, a comparative Monte Carlo study of 
bottle-brush polymers with rigid and relatively long 
backbone lengths {Lf, = 32 to Lf, — 256 monomeric units) 
and flexible side chains of medium length (up to A'^ = 50 
monomeric units) under good solvent and 0-solvent con- 
ditions was performed, using the PERM algorithm. The 
purpose of this study was to investigate the structure of 
such macromolecules and to test physical assumptions 
used in experimental work on related systems to extract 
structural information from scattering data. 

Our main results can be summarized as follows: 

(i) For the chosen side chain lengths, the chosen back- 
bone lengths already are clearly outside of the 
crossover regime from bottle-brush to star polymer 
behavior. Comparing the total scattering function 
Svj{q) of a bottle-brush polymer with and without 
pbcs along the backbone, one does not find any pro- 
nounced effect due to the different conformations 
the chains at the end can assume in the two cases 
(therefore Fig. [9] only shows the scattering for the 
free boundary case). In addition, the range along 
the backbone over which the effect of the proximity 
of the free end of the backbone is felt in the side 
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chain conformations is a few monomer diameters 
only. 

(ii) Corroborating our earlier results [s^l we find scal- 

ing concepts in terms of power laws, blob pictures 
etc. not useful to understand our results. We be- 
lieve that scaling will become useful if the chain 
lengths of the side chains are two orders of magni- 
tude larger; however, this limiting case is beyond 
the reach of either experiment or simulation. 

(iii) Correlations between backbone monomers and 
side chain monomers do not contribute signifi- 
cantly to the scattering , while correlations between 
monomers from side chains anchoring at different 
backbone positions do. As a consequence, the 
standard factorization approximation by which the 
cross-sectional scattering function S^siq) is related 
via Fourier transform to the radial monomer den- 
sity profile p{r), is invalid for most of the momen- 
tum transfer range typically studied. While ex- 
periments typically are done for bottle-brush poly- 
mers with flexible backbones and we deal here with 
the case of rigid backbones only, there is no rea- 
son why approximations that are inaccurate in the 
latter case should become accurate in the flexible 



backbone case, of course. From a detailed anal- 
ysis of the scattering function and radial density 
obtained in the simulation we identify the regime 
where the analysis of the cross sectional scattering 
might be successfully performed to lie at q-values 
smaller than the position of the peak in the curves 
of plots of qSxs{q) vs. q. Here one should fit the 
Fourier transform of a concave form of radial den- 
sity dependence, as given, e.g., by Eq. (fT6|) . 

(iv) It would be desirable to perform neutron scattering 
from bottle brushes where only a small fraction of 
side chains is deuterated. In this way, a more direct 
information on the local conformational structure 
in a bottle-brush could be gained, and more ex- 
tensive comparison with simulations should become 
possible. We also hope that our study will stimu- 
late further experimental work on bottle brushes, 
in particular on the effect of solvent quality. 
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